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In this paper we investigate the spin-current and its shot-noise 
spectrum in a single-molecule quantum dot coupled with a local phonon 
mode. We pay special attention on the effect of phonon on the quan- 
tum transport property. The spin-polarization dependent current is 
Q\ ■ generated by a rotating magnetic filed applied in the quantum dot. 

Our results show the remarkable influence of phonon mode on the 
\ zero-frequency shot noise. The electron-phonon interaction leads to 

sideband peaks which are located exactly on the integer number of 
the phonon frequency and moreover the peak-height is sensitive to 
the electron-phonon coupling. 

g i PACS : 73.23.-b; 85.75.-d ;73.50.Td ;85.65.+h 
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1 INTRODUCTION 

Modern nanotechnology has provided a possibility to fabricate electronic 
devices in which the acting element is a single, organic molecule. This 
leads to a growing interest in the study of transport properties of molec- 
ular devices[TJ 12 Ej. Such a device may be modeled as a 'quantum dot' 
(QD) weakly coupled to the macroscopic charge reservoirs. In addition to 
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their practical applications these artificial , tunable devices such as electronic 
components[T], coulomb blockade structures [I], diodes[5] or switching devices 
with high negative differential resistance , are important for understanding 
the basic physics including the many-body and the size effects. In contrast 
to the semiconductor QDs, molecules ( the linear size is at least one order 
smaller than that of the former) are intrinsically different from semiconductor 
nanostructures. The devices with molecules may lead to new physics espe- 
cially when electrons are added or removed from a single- molecule. Since 
the molecular material possesses much smaller elastic parameters it is very 
easy to excite their internal, vibrational degrees of freedom (phonon modes) 
O 13 IE El EDI UJ] when electrons are incident upon the molecules through 
a tunnel junction. Thus molecules react inevitably back to the tunnel elec- 
trons even at low temperature. This phenomenon has now provoked a large 
amount of experimental [3 |H1 EH HH U21 UH UH US] investigations in the 
problem of transport through mesoscopic systems with electron- phonon 
coupling. Inelastic scattering effects have been observed directly in mea- 
surements of the differential conductance of molecules adsorbed on metal- 
lic substrates with scanning tunneling microscopy (STM)[3]. In a series of 
pioneering experiments by Park et al 8] it was shown that the current 
through a single Cqq molecule was strongly influenced by the vibrational 
mode. Zhitenev et al. [11 have also demonstrated that the low-bias con- 
ductance of molecules is dominated by resonant tunneling through coupled 
electronic and vibrational levels. There has been a number of theoretical 
efforts [T31 HQ HZ1 HH1 [IHl 1201 12H 1221 123] focused on the effects of electron- 
phonon(E-PH) coupling in mesoscopic systems with basic models capturing 
the essential physics and standard methods. Various aspects of the electron- 
phonon/vibron interaction effect on the tunneling through molecule QDs 
have been studied by many authors. 

On the other hand, motivated by the easy control of electron-spin as 
well as the remarkably long coherence time, the spin-polarization depen- 
dent transports in open QDs [201 123 EH] have attracted considerable at- 
tentions. These spin-source devices not only exhibit the new fundamental 
physics but also are of promising applications in the emerging technolo- 
gies of spintronics and quantum information P5 . A pure spin current has 
been reported by direct optical injection without generating a net charge 
current [20.. Theoretically, there are number of mechanisms proposed to pro- 
duce pure spin-current |SU E21 EH] using techniques of ferromagnetic reso- 
nance in a ferromagnetic- normal-metal j31| or electron spin resonance (ESR) 
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in a QD-based system with sizable Zeeman splitting[32j. 

Moreover, a modern trend of the transport studies in mesocopic systems 
is toward not only to considering the transport characteristics of a given 
device but also to examining the noise properties. Due to the discrete na- 
ture of charge carriers, electrical current through a conductor is subject to 
time-dependent fluctuation around its mean value which manifests the con- 
sequence of the quantization of the charge carriers and is usually referred 
to as the shot noise in literature Shot noise defined as the mean-square 
fluctuations of the current flowing through a given terminal at zero temper- 
ature is of great importance and interest because the spectrum of shot noise 
contains additional information about the interactions which the conduction 
electrons undergo [4*0] beyond the mean-current properties and can be used 
to discern different mechanisms resulting in the same mean current. For 
example, shot noise experiments can determine the kinetics of electron and 
reveal information of the correlation of electronic wave functions. Thus it has 
been extensively studied in a wide variety of systems[40j. Furthermore, the 
spin-resolving current correlation is more useful to describe electron correla- 
tion, because the electronic wavepacket with opposite spins is uninfluenced 
by the Pauli exclusion principle and only reflects unambiguous information 
about the interaction. For a two-lead device, correlations can be formu- 
lated by quantities measured at the same lead i.e. the auto-correlation, or by 
quantities measured at the two different leads, namely, the cross-correlation. 
Biittiker 41j pointed out that while cross-correlations of the charge noise can 
either be positive or negative for bosons, they are necessarily negative defi- 
nite for fermions in both the equilibrium and the transport regimes. It has 
been well known [4*2] that anti-bunching in a Fermionic system gives rise to 
negative definite cross-correlation for charge- current. And the conservation 
of charge forces the cross correlation and auto correlation of charge-current 
noise spectra to differ by just a minus sign in a two-probe system. How- 
ever for a pure spin-current the situation is very different because of a lack 
of spin-current conservation due to the spin flipping. The auto-correlation 
of the spin-current is surely positive definite, while the cross-correlation is 
either positive or negative depending on a number of parameter s|34j. The 
zero-frequency shot noise S(0) for a classical conductor|321 is characterized 
by the Poisson value S p (0) = 2e (I) where (/) is the average current, while 
the shot noise in a non-interacting mesoscopic conductor is always reduced 
by the Pauli exclusion in comparison with the Poisson value. Of course, 
shot noise is also influenced by other factors such as electron-electron and 
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electron-phonon interactions. 

Motivated by the achievements in the single-molecule and spin-current 
experiments we in this paper will, based on the ESR mechanism |34[ |S6jof 
generating net spin-current |25j. investigate the electron-phonon effects on the 
spin-current and its shot noise which have not yet been studied theoretically 
for molecular dots. We use the Keldysh nonequilibrium Green function tech- 
nique to calculate the spin current and shot noise through a single-molecule 
coupled to electron reservoirs for the first time and focus on the effect of 
inelastic scattering process. The ESR-type model with a single dispersion- 
less phonon mode is employed to address the vibrational degrees of freedom 
in the molecular dot. All other complexity of real molecular devices, apart 
from interaction with a single longitudinal optical(LO) phonon localized on 
the molecule, is ignored. 



The model system under consideration is illustrated in Fig. 1, which consists 
of a molecule of one relevant level coupled with a single (Einstein) phonon 
mode and two leads which we label as 'left' and 'right'. A time-dependent 
(rotating) magnetic field B(t) is applied in the molecule dot to flip the spin 
of electron. Also a gate electrode is capacitively attached to the dot to tune 
the energy level. The total Hamiltonian of the system is written as 



2 MODEL AND HAMILTONIAN 



H = H L + H R + H P + H D + H'(t) + H T , 



(1) 



where 




(2) 



Hp = w a + a 



(3) 



H D = ^2l(e + crB cos6>) + \(a + a + )]d+d a 



(4) 



(X 



H'(t) = r[exp(— iwt)dfdi + exp(iwt)dtd^\ 



(5) 



with r = Bq sin 9 




(6) 
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The first two terms Hl and Hr of eq.(l) are respectively the Hamiltoni- 
ans for electrons in the left and right non-interacting metallic leads, where 
C^ aa (Ckacr) are the creation (annihilation) operators of electrons with mo- 
mentum k, spin-cr and energy ej~ in the lead a. Here we have set the same 
chemical potential for both leads. The third term H x describes the nondis- 
persive, longitudinal optical (LO) phonon, where w is the frequency of the 
single phonon mode, a + (a) is phonon creation (annihilation) operator. Hp 
and H'(t) correspond to the interaction Hamiltonias between electron and 
phonon in the QD which is subjected to a time-dependent magnetic field 
with uniform strength, 

B(t) = Bo(sin6 cosut, sin 9 sin ut, cos 6 ) 

where B is the constant field strength. Here <i+ (d a ) are the electron-creation 
(annihilation) operators in the QD and e = Eq + eV g is the single energy level 
of the molecule which can be controlled by the gate voltage V g , where e 
denotes the absolute value of electron charge. A is the coupling constant 
between the electron in the molecule dot and the LO phonon mode with 
energy wq . represents the coupling of the molecule with leads, where 
the tunneling matrix elements T^ a transfer electrons through an insulating 
barrier out of the dot. 

2.1 SPIN-CURRENT AND SHOT NOISE FORMULA 

We define the spin-dependent particle current operator in the lead-a as (h = 
1) 

Ja,aa' = ^2 ^ = ^y^^^kagdv' ~ Tka d t C kaa'}, (7) 

k k 

then spin-current operator of spin component a is 

J s a = ^J2 J ^a'Ka' ( 8 ) 

/ 

CTCT 

where cr z is Pauli matrix. 

The spin-dependent current can be computed from current operator Eq. 

(7) 

^At) - (*)> = -£l? te G£^(M) -TLG^,Jt,t)] (9) 
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where the nonequilibrium Green's functions (NEGFs) are defined as 

G^,(t,t) = i(d+(t)C kac ,(t)). 

Using the Keldysh nonequilibrium Green function formalism [37| we obtain 
spin-dependent current 

(10) 

where 

r= E r - 

a=L,R 

is the total tunnel coupling constant which is a function of energy E and 
G r j£; (Ei, E 2 ) are the Fourier transform of the dot-electron retarded (ad- 
vanced) Green's function 

crj>\t, t) = Ti0{±t =f <K (t) , a+(0}> 

in the presence of both the electron-phonon interaction and the tunneling 
coupling between dot and leads, T a (E), the elastic coupling to the a— lead 
(referred to as the line-width function), depends on the hopping strength and 
the density of states p a (E) in the lead-a according to 

T a (E) = 27T Pa (E)\T ka (E)\ 2 . (11) 

In the wide-band limit in which the bandwidth in the leads is much larger 
than both the resonance width and phonon energies, the contact densities 
of states are constant in the region of the resonance. If the hopping matrix 
elements also vary slowly with energy, the couplings with the contacts T a are 
independent of energy either. The Fermi distribution of the lead-a is 

f a (E) = {exp[(£ - v a )/kT] + (12) 

The noise spectra of both charge-current and spin-current can be obtained 
from the correlation S°% between spin-dependent particle currents in lead a 
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and (3: 



naa 
D af3 



A 

J aa (^1/ 



A 

J aa (tl, 



A 

J f3a' (h) 



A 

J /3cr' (£2, 



(13) 



Here (■ • • ) denotes both statistical average and quantum average on the 
nonequilibrium state. 

Briefly, we substitute Eq.(7) and Eq.(9) into Eq.(13) and apply the an- 
alytic continuation so that the exact expression of the zero-frequency, spin- 
dependent correlation can be obtained. It has been shown in Ref. |23j that 
both the cross- and the auto-correlation are needed to characterize the shot 
noise of spin-current for the two-lead system because a spin-current is not 
conserved due to the spin flip induced by the rotating magnetic field. The 
cross correlation shot noise of spin-current is 

Slr = S sp in,i = ((A Jl\— A Jli)(a Jr\— A Jri)) 

= I |f/i(i-/T)r i r ii [|G'^| 2 + |G'^| 2 -2r 2 (|G'^| 4 + |g-^| 4 X]i4) 

and the auto-correlation shot noise is defined by 



Sll = S spini 2 = ((A J LT - A Jli)(A Jl\~ A J Ll )) 

/jrp 
—m - / T ){2[-rir 2 (|qJ 4 + |q T | 4 ) + t l t(\g^\ 2 + |q T | 2 )] 

-r L r R (\G r n \ 2 + \G^\ 2 )}, (15) 

where / T = ft(E) and / ; = f^E-w). 

Once the retarded Green's functions G r aa ;(Ei, E 2 ) are known, the spin 
current and shot noise can be calculated using Eqs. (10) and (13-15). In 
the following we calculate G r ^} (E 1 , E 2 ) with the standard Dyson equation 
approach. To this end we regard the term, which explicitly depends on time t, 
in the Hamiltonian Eq.(l) as the interacting part Hj such that Hq = H — Hi 
. Denoting Green's functions for the Hamiltonian Hn as G 0r ,(e) the full 

° a a v ' 

Green's functions for Hamiltonian Eq.(l) are then calculated from the Dyson 
equation 

G r (E 1 ,E 2 ) = 2nG 0r (E 1 )5(E 1 -E 2 ) 

+ I ^G r (E l , E + E 2 )H'(E)G 0r (E 2 ), (16) 

J 27T 
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where the boldface notation indicates that the electron Green's function in 
the QD and the interacting Hamiltonian H' are the 2x2 matrices in the spin 
space , where the element H' a (E) is the Fourier transformation of H' ai(T2 (t) 
which is seen to be 

H' axC2 {t) = r[e"^4^ + e^d+dt]. (17) 

The full retarded Green's functions of Hamiltonian (1) are then obtained 
from Eq. (16) , after tedious but straightforward algebra, explicitly as[ 



1 - r 2 g[E{) 

G^{E U E 2 ) = 2<k8{E x + aw- E 2 ) % (19) 

1 - r l g[E x ) 

where g{E\) = G® r a (Ei)G%(Ei + aw), a = —a. Using these relations, it is 
straightforward to derive the expression of the spin-current and shot noise 
from Eq.(lO) and Eq. (13-15). 

Since we are interested in the case with strong strength of E-PH in- 
teraction, the Green's function G 0r , (E) can be calculated by performing a 

standard canonical transformation, H = e s He~ s with S = (— )^d+d (T (a + — 

a 

a)[39j. And then all transformed operators are seen to be 



da 


= d a X 


~+ 


= d+ a x+ 




X 


a 






w 



a + — —\^d^d a 
w ^ 



(20) 



(21) 



with 



X = exp[ (a + — a)]. 

Wq 



(22) 
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The electron number operator in the QD is invariant under the transforma- 
tion, 

~H — 

d a d a = d+d a . (23) 
Then the transformed Hamiltonian can be written as 



H = H L + H R + H X + H D + H'(t) + H T (24) 



where 



H L + H R - ^2 £ kC£ aa C kaa 

ktr,a=L,R 

H x = w (a+ - AWi)( fl - AVd+d,) 

a a 

H D = Y][(e + a B cos 6) + \(a + a + - 2— J^d+d^d+d,, 
H'(t) = r[exp(— iwt)d^d± + exp(iwt)d^d^] 

H T = [T ka C+ aa d a X + c.c\. (25) 

k,cr,a=L,R 

The hopping terms between molecule and leads [ Eq. (25)] is modified by a 
factor X, which describes the fact that the electron hopping is accompanied 
by a phonon cloud. Here to avoid unnecessary complication, we consider the 
leads which are unaffected by the phonons. This means that we ignore a 
factor which results from the average of the X operator and does not lead 
to qualitative changes of the tunneling current. The justification for this is 
given in Ref. 16 and Ref. 39. Consequently we have 

Hx ~ Hx 

and 

H = H el + H ph (26) 
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where 

H ei = ^2 £ kCt aa C kaa + ^(e + crB cos 9- A)d+d a 

ka,a=L,R a 

+ yj [TtoC^rf, + c.c] + r[exp(— to£)d+dj. + exp(iwt)d^d^] 

k,rr,a=L,R 

H ph = w a + a 

with A= Due to the E-PH interaction, the single energy level of the 
molecule is renormalized to e' = e— A and then the Green's function ,G 0r , (t) 
, can be decoupled as 

G°JAt) = -i9(t) (Jd a {t)?A0)) (X(t)X + (0)) ph - (%(0)d.{t)) (X+(0)X(t)) p 

(27) 

where 

<l(t) = e^d a e^ a \ 

X(t) = j H pht Xe - iH vht. 

The renormalization factor due to the E-PH interaction is evaluated as |39*] : 

(X(t)X + (0)) ph = e-*« 

(X+(0)X(t)) ph = 

where 

Ht) = g[N ph (l - e iwot ) + (N ph + 1)(1 - e~ lwat )} (28) 

with parameters N ph = exp( ^ o) _ 1 and g = (^) 2 . 

The Fourier transform of the Green's function G°J^, (E) is given by : 

oo 

G%{E) = exp[-g(2N ph + 1)] ^ Ii[2g J N ph (N ph + 1)] exp( to °/ 2 ^ exp"™ " 

£=— oo 

*[(1 - (v»G w (£ - too) + K ff > G aa {E + too)] (29) 

where G aa (E) is the retarded Green function corresponding to the time- 
independent part of the new Hamiltonian H e i , the index I indicates the 
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number of phonons involved, and (n^) is the time- averaged electron occu- 
pation number in the molecule. Here we consider the case of an "empty" QD 
i.e. (rid,o) =0. With a little algebra we find 

G ™' {E) = E-(e + aB Ze- A) - (3 ° } 

where the retarded self-energy due to the tunneling into the electrical leads 
are given by 

In the wide-band limit, the level shift A.(E) can be neglected and the linewidths 
are energy independent constants. Thus the retarded self-energy can be ex- 
pressed as 

KAE) = ~ l lf- (32) 

The Fourier transform of the full Green's function given by Eq. (27) can be 
obtained as 

G°JAE) = exp[-g(2N ph + l)} (33) 

pWoW/2 



*J2l l [2g^N ph (N ph +l)]- 



E-(e + aB cos 9- A) - lw + \T 

l = — OO v 'A 



3 NUMERICAL RESULTS AND DISCUS- 
SIONS 

We now present the numerical results of the spin current and zero-frequency 
shot noise. For simplicity, we consider the symmetric tunnel-coupling be- 
tween the molecule and the two leads i.e. = = ^ and further assume 
that the energy level of the molecule which is controlled by the gate volt- 
age V g such that e(V g ) = eo + eV g , where eo denotes single-electron energy 
in the molecule in the absence of the gate voltage . The phonon energy is 
chosen as the energy unit throughout the rest of the paper. We also set 
h — e — 1. In fig. (2) we plot the spin current I s versus the parameter r 
(Fig. (2a)), the gate voltage V g (Fig. (2b)) and frequency w of the rotating 
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magnetic field respectively (Fig. (2c)) at zero temperature T = 0. The pa- 
rameter values used in Fig. (2a) are such that T = 0.04, V g — 0, g — 0.5 2 
and 9 = 88°. For comparison, we also plot the spin current in the absence 
of E-PH interaction(dashed line). It is clearly shown in Fig. (2a) that the 
E-PH interaction results in the shift of resonance peaks. Fig. (2b) shows 
the spin-current I s as a function of the gate voltage V g for different r with 
T = 0.04, (7 = 0.6 2 ,w = O.land 6 = 88° . In the presence of E-PH coupling 
the overall spectrum is shifted by a quantity A= ^- toward the negative 
gate voltage region. In addition to the main peak which is related to the 
molecule energy level, small satellite resonant peaks appear at the positive 
energy side. At zero temperature only phonon-emission processes are allowed 
since before tunneling the phonon state is vacuum, which explains why the 
satellite peaks are located at the positive energy region. Moreover, the height 
of the satellite peaks is much smaller than the main resonant peak because 
of the suppression by the E-PH coupling. Fig. (2c) display the spin current 
I s versus frequency w of the rotating magnetic field for various coupling con- 
stants g = 0, 0.3 2 , 0.5 2 , 0.7 2 (correspondingly the solid, dotted , dashed and 
the dash-dotted curves),with T = 0.04, V g = 0, r = 0.08,0 = 88°. In the 
absence of the phonon mode, there exists only one resonant peak. The E- 
PH coupling typically leads to the new satellite resonant peaks. It is also 
seen that the heights of satellite peaks increase with the coupling constant 
g. The positions of the side peaks are located in w = nw {n = 1, 2, 3, 4 • • • ) . 
Finally, the plot of spin current I s as a function of E-PH coupling constant 
g ( Fig. (2d)) shows a double-maximum. Fig. (3) shows the dependence of 
cross-correlation shot noise on the parameter r (Fig. (3a)), the gate voltage 
V g (Fig. (3b)) and frequency w of the rotating magnetic field(Fig.(3c)) respec- 
tively with the same parameter values as in Fig. (2). We see from Fig. (3) that 
the cross-correlation shot noise of the spin-current displays very different be- 
havior from the spin-current itself, therefore the measurements of the shot 
noise spectrum can provide more information of the transport properties in 
mesoscopic systems. From Fig. (3a) it is observed that in the absence of E-PH 
interaction the cross-correlation shot noise can be either positive or negative 
as r is changed due to the competition [27j between the cross correlations of 
electrons with parallel and antiparallel spins. The E-PH interaction has con- 
siderable influence on the cross-correlation shot noise that the negative shot 
noise is substantially suppressed and, even vanishes. Fig. (3b) shows that 
with r increasing, the cross- correlation shot-noise spectra exhibits two extra 
peaks located symmetrically around the position of the main peak, and in 
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the case of stronger magnetic field, the cross-correlation shot noise turns to 
negative for some gate voltages while it remains positive for the other gate 
voltages. From Fig. (3c) the oscillating behavior of the shot noise between 
positive and negative values is observed and is due to the photon assisted 
process. We finally plot the auto-correlation shot noise versus the gate volt- 
age V g ( Fig. (4a) ) and frequency w ( Fig. (4b)) in Fig. 4. The auto correlation 
shot noise spectra are also shifted and satellite peak appears compared with 
the absence of the E-PH interaction. Moreover, the auto-correlation of the 
spin-current is surely positive definite which is different from the situation 
for cross-correlation. We see that the shot noises possess reach informations 
of the phonon effect on the quantum transport through molecular devices. 

4 CONCLUSIONS 

In summary we have shown for the first time that the shot noise of spin 
current in a single-molecule quantum dot can be significantly affected by the 
phonon mode and can provide a beneficial information to improve the un- 
derstanding of transport properties through the molecular QDs. It is shown 
that in addition to the shift of the resonant-peak position associated with 
the level of the dot, satellite peaks emerge at integer number of the phonon 
frequency. The E-PH coupling even can reverse the sign of the zero-frequency 
cross shot noise. 
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Figure Captions: 

FIG.l. Schematic diagram of a molecule quantum-dot system. 

FIG. (2a)The spin current I s versus the parameter r ( T = 0.04, V g = 
0,g = 0.5 2 , 9 = 88° with frequencies of magnetic field w = 0.05,0.25,0.7 
and 1.5 ) for the two cases with (solid line) and without the phonon mode ( 
dashed line) as comparison. 

FIG. (2b) The spin current I s as a function of the gate voltage V g ( Y — 
0.04,# = 0.6 2 ,w = 0.1 and 9 = 88° with the parameters r = 0.02,0.04,0.06 
and 0.08) for the two cases with (solid line) and without the phonon mode ( 
dashed line) as comparison. 

FIG. (2c) The spin current ^ versus frequency w of the rotating magnetic 
field with T = 0.04, V g — 0, r — 0.08, 9 = 88° for various coupling constants 
g = (solid curve) g = 0.3 2 (dotted curve) g = 0.5 2 (dashed curve) g = 0.7 2 
(dash-dotted curve). 
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FIG. (2d) The spin current I s as a function of coupling constant g . 

FIG. (3a) The cross- correlation shot noise as a function of the parameter 
r (solid line: with phonon mode; dashed line: without phonon mode). 

FIG. (3b) The cross-correlation shot noise as a function of the gate voltage 
V g (solid line: with phonon mode; dashed line: without phonon mode). 

FIG. (3c) The cross-correlation shot noise versus frequency w for 
different coupling constants: g — (solid curve) g = 0.3 2 (dotted curve) 
g = 0.5 2 (dashed curve) g = 0.7 2 (dash-dotted curve) with T = 0.04, V g = 0, 
r = 0.08,0 = 88°. 

FIG. (4a) The auto-correlation shot noise versus the gate voltage V g with 
r = 0.02, 0.06, O.land 0.2 (solid line: with phonon mode; dashed line: without 
phonon mode). 

FIG. (4b) The auto-correlation shot noise versus frequency w with the 
same parameters as in Fig. (2c) for various coupling constants g = (solid 
curve) g = 0.3 2 (dotted curve) g = 0.5 2 (dashed curve) g = 0.7 2 (dash-dotted 
curve). 
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Fig.(l) 
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